clear;
% 读取图像
I = imread('fig6.tif');

% 添加椒盐噪声
J = imnoise(I,'salt & pepper',0.02);

% 使用 3x3 平均滤波器
K1 = filter2(fspecial('average',3),J);

% 使用 5x5 平均滤波
K2 = filter2(fspecial('average',5),J);

% 使用 3x3 中值滤波器
K3 = medfilt2(J,[3 3]);

% 使用 5x5 中值滤波器
K4 = medfilt2(J,[5 5]);

% 显示结果
subplot(2,3,1),imshow(I),title('原图');
subplot(2,3,2),imshow(J),title('添加椒盐噪声');
subplot(2,3,3),imshow(uint8(K1)),title('3x3 平均滤波器');
subplot(2,3,4),imshow(uint8(K2)),title('5x5 平均滤波器');
subplot(2,3,5),imshow(K3),title('3x3 中值滤波器');
subplot(2,3,6),imshow(K4),title('5x5 中值滤波器');